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Ph 1 INTRODUCTION 



By considering that the distribution of measured light along any galactocentric 
radius of an elliptical galaxy can be well represented by a i?^/" Sersic law, we propose 
in this paper a "discrete ordenate" method, which, for any value of n, allows an 
explicit expression for the luminous density, p{r), to be found that can be numerically 
evaluated to any required precision. Once we have this semi-analytical expression for 
the spatial density, p{r), the mass distribution, M{r), the potential, (/'(r), and the 
velocity dispersions, cr^(r), in the space and in the observational plane, (7p{R), can be 
computed in a straightforward manner. 



Key v^rords: galaxies:elliptical and lenticular, cD 
kinematics and dynamics. 
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Analysis of the distribution of the light intensity over the 
surface of an elliptical galaxy is fundamental to determining 
its internal structure. From the beginning it was pointed 
out that this intensity distribution over a galactocentric 
radius — for simplicity along the major axis of the observed 
ellipse — was similar for all galaxies, as well as for the bulges 
of spiral galaxies. This distribution was firstly represented 
by the Hubble-Reynolds law I{R) oc (r + ro)~^ (Reynolds 
1913; Hubble 1930), and some time later by the "universal" 
de Vaucouleurs law. However, later observations (Davies et 
al. 1988; Caon, Capaccioli & D'Onofrio 1993; Andreakis, 
Peletier & Balcells 1995; Graham et al. (1996); Binggeh & 
Jergen 1998) did show some discrepancies with the R^''^ law, 
and led to a more general form with B}'" , generally with 
n > 1 (Sersic 1968). This bi-parametric expression accounts 
much better for the observed intensity profiles of elliptical 
galaxies and bulges of spiral galaxies. 

In this paper we present a careful study to determine 
the 3D spatial distribution of emitting sources from any ob- 
served _R^' " Sersic profile, and via a mass-to-luminosity re- 
lation we derive the corresponding dynamical properties. 



2 THE SERSIC PROFILE 

For an elliptical galaxy the "light profile" which describes 
the light intensity distribution over any galactocentric radius 
R{oi) corresponding to the galactic longitude a can be well 
approached by the Sersic profile. 



I{R{a)) = /(0)exp[-(i?(a)/i?o(Q))" 



(1) 



where the value of n is the same for all the galactocentric 
radius R{a). The scale length, Ro{a), will have different 
values for different a according with the elliptical shape of 
the galaxy. If a and b are, respectively, the major and minor 
semi-axis of a given isophote (e = 1 — b/a is the ellipticity) 
and Ro = Ro{0) is the scale length of the light profile along 
the major semi-axis, we then have. 



Ro{ct) = Ro- 



(1-^) 



'a-h(l-e)2 



(2) 
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Equations (1) and (2), with n independent of a, indicate 
that all the isophotes (curves with the same value of the light 
intensity, I) are ellipses with the same ellipticity, e. 

In general, the observational intensity profiles are mea- 
sured along the major axis of the galaxy, where a = and 
therefore we note I{R) = I{R{0)). However, other parame- 
ters have been proposed for R (radius) for identifying each 
observed isophote and measuring the intensity profile. For 
example, if at and bk are the major and minor semi-axes of 
each isophote, respectively, instead of R — ajc we can use 
the value of b^, the mean value (aj, + bfc)/2 or the geometric 
mean (ytifefcfc). However, each of these options leads to the 
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same conclusions. We might have only small differences from 
the practical point of view. Therefore, purely for sake of sim- 
plicity we carry on with our work using the semi-major axis 
as the "radius parameter", R, to measure the distribution 
of luminous intensity over the galaxy. 

The total energy or luminosity enclosed by an isophote 
defined by the value R of its major axis is given by. 



L{R) 



1 



-R{«) 



da / R'{a)I{R'{a))AR!(a), 



(3) 



R[a) being the value of the radius vector of the points of the 
ellipse-isophote defined by R that corresponds to the galac- 
tocentric longitude a. If the brightness distribution follows 
the Sersic law. Equation (1), and taking into account that 
for the isophote defined by R we have R{a)/ Ro{a) — R/Ro, 
then Equation (3) reduces to 



L{R) = I{0)RItt{1 - e)2n7(2n, {R/Ro)-), 



(4) 



7(2n, {R/RoY'") being the incomplete gamma function cor- 
responding to the gamma function r(2n) for the value 

(i?/i?o)^/". 

The total observed luminosity is given by Equation (4) 
when R -^ oo: 



7(0)K7r(l-e)2nr(2n). 



(5) 



But instead of using the natural scale length Ro to de- 
scribe a profile I[R) measured along the semi-major axis of 
the isophotes, the value Re of R, for which the correspond- 
ing isophote encloses the half of the total luminosity is cus- 
tomarily used. Obviously, the other half corresponds to the 
luminosity over the isophotes with R > Re- This definition 
of Re can be written as 



-r{2n)=^{2n,{Re/Ro)-). 



(6) 



This equation gives a value of Re/ Ro for each n. If we 
define k — {Re/RoY , the values of k that satisfy Equation 
(6) can be easily computed for any value of n. For n > 1, 
they can be obtained (with an error smaller than 0.1%) by 
the relation fc = 2n - 0.324 (see Ciotti 1991). 

The use of the variable k and the characteristic radius 
Re allow us to write the Sersic profile. Equation (1), in its 
usual form (see Fig. 1): 



/(i?) = /(0)exp[-fc(i?/i?e 



(7) 



Likewise, the luminosity distribution in Equation (4) 
becomes. 



L{R) = I{Q)RlA^ - e)^l{2n, (R/R 
and the total luminosity is, 
L = I{0)Rln{l~e)^r{2n). 



■ej ), 



(8) 



(9) 



In Fig. 1 we plot, together with the observed intensity 
profile I{R), the normalized observed luminosity, L{R)/L, 
for different values of n. 



3 THE ABEL INVERSION FOR THE SERSIC 
PROFILES 

The problem that we want to solve is to determine the 
galactic density of luminous sources which can explain the 
observed intensity profile, I{R), given by Equations (1) or 
(7). This intensity at a given position is originated for by 
the emission of all the stars that are located along the line 
of sight. For the case in which the emitting sources of the 
galaxy have a spherical distribution of density, Pl(t), the 
observed isophotes will be circular and the corresponding 
intensity will be given by 



,/„x I , s 2rdr 



(10) 



where R is the radius of each observed isophote. The prob- 
lem that we have to solve is to determine the density of lu- 
minosity, pL{r), for a known intensity profile, I{R)- In math- 
ematical terms the problem consists in the inversion of the 
Abell integral transform with an infinite interval (Equation 
[10]). 

If the distribution of emitting sources is not spherical, as 
happens in the case of elliptical galaxies, the problem is not 
so simple. Here we can admit that the density of sources 
pL is constant over concentric and coaxial ellipsoidal sur- 
faces which can be bi- or triaxial. Let d be the parameter 
whose value characterizes each of the ellipsoidal shells in 
the space. It can be the value of one of its three semi-axis, 
a mean (linear or quadratic) value of these axes, or any 
other parameter with linear dimensions. Then the density 
of luminous sources, pL^d), will be a function of d. If these 
ellipsoidal shells are homologous (i.e. if the ratio between 
their semi-axes are the same for all of them) the correspond- 
ing isophotes over any observation plane are concentric and 
coaxial homologous ellipses (all have the same ellipticity) . 
Taking D as a parameter to describe the observed isophotes 
[D can be the semi-major axis or any other linear dimension 
of these isophotes as we discuss in Section 2), we can write 
an equation similar to Equation (10) (see Lindblad 1956), 
which relates the intensity, I{D), of each isophote with the 
density of emitting sources, pL{d). Now d plays the same role 
as r in the case of a spherical distribution of sources, but d 
is no longer a vector position in the space. The choice of the 
parameter D to characterize each of the observed isophotes 
and to describe the distribution of the intensity, I{D), over 
these depends on the choice of the parameter d to describe 
the structure and emission density, pL (d) , of each ellipsoidal 
shell that forms the galaxy. To find the relation between 
D and d we need a detailed geometrical description of the 
shape, position and projection of these ellipsoidal surfaces 
(see Simonneau, Varela & Munoz-Tunon 1998). 

However the main goal of this work is to describe in a 
semi-analytical way, some mathematical properties of a den- 
sity pL{r) obtained from an inversion of the Abel integral 
transform (Equation [10]) of an observed radial profile, I{R), 
which follows the Sersic law (Equation [7]). These mathe- 
matical properties are the same for the case of a spherical 
galaxy where the emitting shells are defined by r and the 
circular isophotes are defined by its radius, R, which, for 
the case of an ellipsoidal galaxy where the emitting shells 
are defined by the parameter d and the elliptical isophotes, 
are defined by another parameter, D. In this paper, then. 
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we shall consider only the case of the spherical galaxies by 
means of Equation (10). But we would like to emphasize 
that many of the conclusions presented in this paper can be 
applied, with minor changes, to the case of elliptical galax- 
ies, once we find out how to define the d and D parameters 
correctly. 

In terms of the new variables z = (R/Re) and s = 
(r/Re), the theoretical inversion of the Abel integral trans- 
form (Tricomi 1985), equation (10), leads to 



RepL^s) = 



dl{z) 



dz 



dz Vi5- 



(11) 



where I{z) is given by the Equation (7). It will then be 



Pl{s) 



1 k J(0) 

TT n Re 



or otherwise 

, , fc/(0) i_i 



exp[— fca:"] z' 



exp[— fcs"i] 



1 



-dz, 



y? 



dt. 



(12) 



(13) 



The argument in the integral of Equation (13) is sin- 
gular at t = 1, although it is integrable for any value of n. 
Only when s = does the integral diverge for n < 1 be- 
cause of the infinity in the upper limit. But for these values, 
because of the multiplicative factor s^'"~^, Pl{s) becomes 
zero for s = 0. When s increases, Pl{s) also increases until 
some critical value and then decreases monotonically. This 
means that some layers are more dense than the internal 
ones, and that consequently the system is gravitationally 
unstable. Hence, models with n < 1 are not adequate for 
representing stable galaxies. 

When n — 1 the intensity profile, I{R), takes an expo- 
nential law and the integral in Equation (13) is the Ko{x) 
modified Bessel function, which has a logarithmic discon- 
tinuity for X = (see Abramowitz & Stegun 1968, pages 
374-376). Hence, for n = 1 the spatial density becomes: 



n Re 



(14) 



For n > 1 the integration of Equation (13) is not possi- 
ble analytically, except for s = 0, where 



dt 



Vi^" - 1 



^2^^^r- 



2n 



-) 



(15) 



and B{z,'w) is the complete beta (Euler) function (sec 
Abramowitz & Stegun 1968, page 258). 

However, the form of the integral in Equation (13) has 
been used by many authors to obtain selected approxima- 
tions for Pl{s). First, for the de Vaucouleurs case n = 4, 
by Poveda, Iturriaga & Orozco (1960), Mellier & Mathez 
(1987), and for the Sersic profiles (Gerbal et al. 1997). Also, 
numerical computations of Pl{s) have been carried out by 
Young (1976) for the de Vaucouleurs case and by Ciotti 
(1991) and Graham & CoUes (1997) for the Sersic cases, 
respectively. In this article we propose an analytical expres- 
sion for pl (s) that allows an easy computation of the mass 
and gravitational potential to any required precision. 

As the argument of the integral in Equation (13) can be 
integrated for any value of s, even considering the singularity 
for f = 1, it seems possible to write this integral with a 
new variable, x, such as its argument does not show any 



discontinuity. We propose t — 1/(1 
(13) for the density becomes 



„2U/(n-l) 



Equation 



Pl(s) = - 



TV Re n — 1 



cxp[— fcS" (1 — x^) 



(1-^2)^^ 



a;da;.(16) 



We can then perform this integral formally by means of 
a Gaussian numerical integration to get 



Pl{s) = 



where 



k /(O) 2 



TV Re n 






J=l 



(1-^?)^ 



Pj = Wj 



l-{l-xf)- 



(17) 

(18) 
(19) 



{j — 1,2, ..., Nap) are the abcissae and the 
integration weights, respectively, in the interval (0, 1). Nap 
is the order of approximation, i.e. number of abscissae and 
weights needed to compute the integral in Equation (13) 



numerically. The standard tables of abcissae, 



and in- 



tegration weights, Wj, correspond to the interval (—1,-1-1) 
(Abramowitz & Stegun 1968, pages 916-919). Therefore, the 
relation of the values of Xj and Wj for the interval (0,1) are 
Xj — {l + Xj)/2 and Wj=Wj/2. In summary, by selecting dif- 
ferent sets of abcissae, Xj, and integration weights, Wj, we 
can compute Xj and the coefficients pj by Equations (18) 
and (19), respectively, and determine to any required pre- 
cision the luminous density Pl{s) from Equation (17). In 
Table 1, we have listed the abscissae and weights for the 
Nap = 1,2 and 5 cases that we discuss below. 

In the first place, we must point out that there are few 
characteristic quantities in relation to the distribution Pl{s) 
that can be determined theoretically, i.e. independently of 
the numerical integrations. One of these is the integral in 
Equation (15), which determines the asymptotical behaviour 
of Pl{s) when s -^ 0. From Equations (15) to (17) it is 
necessary that 



4n 



Nap 



Pj ~ B(t 



2n 



(20) 



where the coefficients pj are those given by Equation (19). 
The other characteristic quantities are the different moments 
of the density distribution, pL(r), i.e. 



pL(r)r^dr. 



(21) 



The first one (m = 0) loads to the synthetic value of 
/(O). To satisfy Equation (10) for R — with a density law 
as in Equation (17) it is necessary that 



A, ~ 2 



n n 
2 



2n 



(22) 



j=i 



The second one (m = 1) leads to the value of the grav- 
itational potential at the centre of the system {r — 0). This 
value, (jf'(O), is known a priori (see bellow Equation [31]). To 
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Table 1. Abscissae, Xj, and weights, Wj, for the Gaussian 
integration in the interval (0,1) for three different approxima- 
tions, Afap=l,2 and 5. 



a-p 


= 1 


0.5 


1. 


'ap 


= 2 


0.211325 


0.5 






0.788675 


0.5 


ap 


= 5 


0.046910 


0.118464 






0.230765 


0.239314 






0.5 


0.284444 






0.769235 


0.239314 






0.953090 


0.118464 



reach this theoretical value with a density law, as in Equa- 
tion (17), is necessary that 






n-1 
2n 



(23) 



The third one (m = 2) leads to the value of total lumi- 
nosity, L. To have the same theoretical value as in Equation 
(9) it is necessary that 



E Pj 



4 2n 



(24) 



directly from Equations (20), (22), (23) and (24) in such a 
way that, even if between these first orders the computed 
Pl{s) and the correct solution there are any differences, we 
are sure that the most important parameters of the problem 
are computed exactly. 

For most applications where great accuracy is not 
needed, the approximation Nap ~ 5 can be sufficient. In 
case of much higher accuracy the Nap = 10 and Nap = 20 
approximations can satisfy any requirement. 

In Fig. 2, we plot Pl{s) as a function of s = r/R^ for 
different values of n. 



Thus, to test as a first step the accuracy of our ap- 
proximation for Pl{s) in Equation (17), we have computed 
the left-hand side of Equations (20), (22), (23) and (24) nu- 
merically; and in Tables 2, 3, 4 and 5 we have compared 
these with the theoretical values of their right-hand side to- 
gether for the different approximations Nap~l, 2, and 5, 
and for different values of n (2, 3, 4,.. .,10). As can be seen, 
the approximation Nap = 5 leads to results of high quality. 
According to the results we have that for Nap ~ 10,20 and 
40, the variation is insignificant with respect to the results 
for Nap = 5. 

We now discuss the accuracy of each approximation 
(A'ap=l,2,5,10,20,40) over the entire radial density distri- 
bution, p{s). As a reference we take the density p4o(s) cal- 
culated with Nap ~ 40. The relative variation between this 
and p2o{s), calculated with the approximation Nap = 20 is 
always smaller than 0.01%, i.e. negligible to the first four 
significant figures. This difference increases by one order of 
magnitude for the Nap = 10 approximation. The relative 
difference between pio(s) and p4o(s) is always smaller than 
0.1%. For the case of Nap ~ 5 the difference between P5(s) 
and p4o(s) never reaches 1%. For the case of Nap = 2 we can 
get differences of around 10% only for low values of n (less 
that 5), and for large values of s (greater than 10). This 



relative difference increases significantly for Na 



1. So, 



if we consider as not significant the difference between the 



Nap = 40 and iV„ 



20 cases, we can take as the errors of 



each approximation the differences given above. In general. 
Nap = 2 is a good approximation for getting a complete 
first-order description. However, concerning the quality of 
the first approximations (A^ap = l,2) we cannot consider the 
Gaussian numerical integration as the best method for ob- 
taining the values of pj,Xj. These values must be obtained 



4 DYNAMICAL PROPERTIES FOR THE 
SERSIC PROFILE 

Once we have the density of luminous sources, pl{s), we 
can study some other quantities in relation to the dynamical 
state of the galaxy. 



4.1 The mass distribution 

The density Pl{s) in Equation (17) refers to the density of 
luminous sources. In cases in which the mass-to-luminosity 
ratio, T, is the same in all the points of the galaxy, the 
total mass, M, will be given by Equation (9) multiplied by 
this ratio T. Then, we can use p{s) for the mass density 
distribution with the form 



p{s) = PoP{s), 

where (with M = Ti), 

L 2 A;2"+i 



po^T 



and 



P{s) 



7r2i?i n - 1 r(2n -f 1) 



;r^y^Pjexp[-Ajfcs~]. 



(25) 



(26) 



(27) 



The corresponding mass distribution, M{s), 

M{s) =Ri I 47rs'^p(s')ds, (28) 

Jo 

is given by 
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Table 2. Asymptotic coefficient of Pl(s) for s ^ 0, Equation 
(20). We show the exact value of B(i, "g" ) ^^^d ^^^ corre- 
sponding computed values " V^ .J^'' pj, Equation (20), with 



Na 



=1,2,5. 



n 


Bih^) 


Nap = l 


Nap =2 


Nap =5 


2 


5.2441151 


4.837945 


5.255770 


5.244083 


3 


4.2065460 


3.945576 


4.204081 


4.206550 


4 


3.8558066 


3.643523 


3.850307 


3.855832 


5 


3.6790940 


3.490923 


3.672538 


3.679124 


6 


3.5725536 


3.398728 


3.565529 


3.572583 


7 


3.5012896 


3.336962 


3.494024 


3.501318 


8 


3.4502622 


3.292684 


3.442860 


3.450289 


9 


3.4119198 


3.259382 


3.404436 


3.411943 


10 


3.3820539 


3.233423 


3.374518 


3.382076 



Table 3. Synthesis of the observed central intensity /(O). We 
show the theoretical value of ^ 2^ together with the cor- 

Z An <=' 

responding numerical computation of Y^ .^^ pj /\j, Equation 
(22), with Afap=l,2,5. 



7r n— 1 
2 2n 



Nap 



1 



Nap 



Nap ■■ 



2 
3 

4 
5 
6 
7 
8 
9 
10 



0.3926991 
0.5235988 
0.5890502 
0.6283185 
0.6544985 
0.6731984 
0.6872234 
0.6981317 
0.7068726 



0.453557 
0.569495 
0.620693 
0.649734 
0.668478 
0.681587 
0.691273 
0.698723 
0.704633 



0.397601 
0.537842 
0.603678 
0.641712 
0.666468 
0.683865 
0.696759 
0.706698 
0.714593 



0.392702 
0.524518 
0.590514 
0.629974 
0.656182 
0.674842 
0.688799 
0.699631 
0.708280 



M{s) 



Nap 



Pj 



M 7r(n - l)r(2n) -:1^ A^' 



'—-f{2n + l,Xjks-), (29) 



where 7(2n + 1, Ajfcs") is the incomplete gamma function 
corresponding to the gamma function r(2n-|-l) for the value 
\jks~ . Because the Gaussian values of \j and pj satisfy 
Equation (24) very accurately we are sure about the correct 
normalization in the M{s) distribution. 

The normalized mass distribution, M{s)/M, is shown in 
Fig. 2 for the different values of n. For n — 1 the distribution 
of mass has been computed directly from the corresponding 
density given by Equation (14). 



4.2 The gravitational potential 

For a spherical system with density distribution p{s) (Equa- 
tion [25]), the gravitational potential is given by 

(t)(s) ^ -AttGRIpo{- p(s')s'^ds'+ / p{s')s''^ds')l,{30) 

where, for the Sersic law, po and p{s) are given by Equations 
(26) and (27), respectively. 

At the centre, for s = 0, we have 



^(o) = -(^)^fc"^^"^ 



Re 'it r(2n)' 
© 0000 RAS, MNRAS 000, 000-000 



(31) 



which is the correct value known a priori (see Ciotti 1991). 
For the potential we then have 



(s) 



2n(j!>(0) 



(n - l)r(n + 1) ^ A, 



^ ;iTr'>("+l.^i^s")~(-^)-— — ,(32) 



Re ' s M 



where M{a)/M is given by Equation (29). As before, the 
accuracy of the sum X^ .J!f Pj/A"^^ in Equation (23) leads 
to the correct value for 4>{Q) in Equation (32). 

Curves for the potential (^{s) for different values of n 
(n = 1, 2, ..., 10) are shown in Figure 3. 



4.3 The total potential energy 

With the above expression for p{a), M{s) and 0(s) it is very 
simple to compute the total potential energy. 



to get 
W = T 



Anr p{r)(f>{r)Ar, 



2 (jL 2 

-rT""- 



(33) 



(34) 



The values of the parameter w'^ , computed numerically 
from Equation (33), are given in Table 6 for the different 
values of n. 
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Table 4. Synthesis of the theoretical gravitational potential 
at r = 0. We show the theoretical value of "2" together with 
the corresponding numerical computation of V^ .^^ pj/\ , 
Equation (23), with Nap=l,2,5. 



n 


n-l 
2n 


Nap =1 


Nap =2 


Nap = 5 


2 


0.25 


0.255126 


0.246949 


0.249997 


3 


0.3333333 


0.369898 


0.327389 


0.333356 


4 


0.375 


0.422952 


0.370122 


0.374965 


5 


0.4 


0.453484 


0.396269 


0.399942 


6 


0.4166667 


0.473327 


0.413849 


0.416597 


7 


0.4285714 


0.487259 


0.426461 


0.428498 


8 


0.4375 


0.497580 


0.435943 


0.437496 


9 


0.4444444 


0.505533 


0.443330 


0.444371 


10 


0.45 


0.511849 


0.449245 


0.449928 



Table 5. Synthesis of the theoretical value of the total mass. 

together with the cor- 



We show the theoretical value of ^ ^ — 
responding numerical computation of Y^ 
tion (24), with Nap = 1,2,5. 



Nap ^,/x2n4 

1 pj/^j 



Equa- 



n 


7r n— 1 
4 2n 


Nap = l 


Nap =2 


Nap = 5 


2 


0.1963495 


0.143508 


0.208820 


0.196357 


3 


0.2617994 


0.240256 


0.265076 


0.261793 


4 


0.2945431 


0.288208 


0.294168 


0.294525 


5 


0.3141592 


0.316511 


0.312086 


0.314166 


6 


0.3272492 


0.335146 


0.324230 


0.327259 


7 


0.3365992 


0.348335 


0.333001 


0.336611 


8 


0.3436117 


0.358159 


0.339631 


0.343623 


9 


0.3490659 


0.365757 


0.344818 


0.349077 


10 


0.3534292 


0.371810 


0.348987 


0.353441 



The ratio Re/w^ can be considered as a mean grav- 
itational radius, re, such that W = T^GL^ /ra- That is, 
in some way tq = Re/w^ accounts for the mass concen- 
tration degree in the inner parts of the system. But from 
another point of view we can interpret {TGL/Re)'w^ as a 
quadratic mean velocity such that the corresponding kinetic 
energy, T = (l/2)M(TGL/i?e)w^, satisfies the virial theo- 
rem, W = 2T. We should now examine the possibility of 
measuring this quadratic mean velocity. 

In stationary spherical galaxies we can assume that 
there are no organized motions of the stellar populations 
as a whole, i.e. rotation, expansion or contraction; hence at 
each point in the galaxy we will have a quadratic mean ve- 
locity corresponding to the local velocity dispersion { crs^(r) 
). Indeed this one is a double mean. First, it is a local 
quadratic mean of the three components of the velocity dis- 
persion (Ts{r)= {ar{r),a^{r),ae{r)). However, to be able to 
apply the virial theorem we need the necessary condition of 
quasi-isotropy in the velocity dispersion tensor. In practice 
we can admit isotropy, i.e. ar{r) = cr^(r) — (Jg{r) = (Js{r), 
to be able to apply the virial theorem with a quadratic mean 
velocity: 



{'Ts\r))=3 



/y°°crs(r)r^p(r)dr 



3a, 



(35) 



If observations of this quadratic mean velocity (disper- 
sion velocity) are available, and if we admit a "Sersic model" 



for the galaxy and consequently compute, for any observed 
I{R) profile, the potential energy (Equations [34]; see Table 
6 for the factor w^ corresponding to each n), we could write 
the virial theorem in the form 



M3ai 



2(GL. 2 



(36) 



where the total luminosity, L, is another measured quan- 
tity. Therefore, we can estimate the value of the mass-to- 
luminosity ratio, T. This method for the determination of 
the masses of elliptical galaxies was proposed, for the case 
n = 4 (the de Vaucouleurs law) by Poveda (1958). 

However, we do not have access to measurements of 
the velocity dispersion, (Js{r), in the space. We can measure 
only the velocity dispersion, Gp{R), in the observation plane. 
Theoretically, these measurements can be represented by 



l(R)<Tp{R) 



{r)pL{r) 



2r 



\pr^ 



i?2 



dr, 



(37) 



which corresponds to the integral along the line of sight 
of the spatial velocity dispersion, ai,(r)^ weighted with the 
density. The quadratic mean value of ap{R) is 

j;^ 4{R)dL{R) _ , 



i^p) 



(38) 



From Equations (35) to (37) and from the definition 
of I{R) in Equation (10), the condition that a^ — a'j is 
satisfied. This gains us access to the quadratic mean spatial 
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Table 6. The potential energy, w'^, measured by taking 
(T^GL^/Re) as unity in Equation (34) for different values 
of n. 



1 


0.31426 


2 


0.30973 


3 


0.31856 


4 


0.33615 


5 


0.35983 


6 


0.38897 


7 


0.42349 


8 


0.46363 


9 


0.50981 


10 


0.56260 



velocity dispersion, 3a^, that in the virial theorem leads to 
the mass-to-luminosity ratio, T. 



4.4 The observed velocity dispersion 

The measurement of the velocity dispersion on the obser- 
vation plane, Equation (38), needs an integration over the 
entire galaxy — mathematically, over the entire radial inter- 
val, i.e. (0, (X)). In general, this condition cannot be satis- 
fied and we can measure only I{R)a^{R) in a radial interval 
(0,7?m) which does not contain the entire luminosity, L, of 
the galaxy. This can lead to an incorrect value of a^ and 
therefore a wrong estimate of the mass-to-luminosity ratio, 
T, from the virial theorem. 

However, this lack of spatial coverage in the observa- 
tions can be compensated thanks to our galaxy model. Once 
we have the density distribution law. Equation (17), and the 
potential distribution. Equation (32), the spatial velocity 
dispersion, ^^(r) (assuming isotropy in the velocity distri- 
bution function), satisfies the Maxwell-Jeans equation. 



— p(r)3o-3(r) = -p(r) — 0(r), 



(39) 



and, with the natural boundary condition p{r)a1{r) — > for 
r -^ oo, we have 



d 



p(r)3cr,(r) = / p{r')-—d){r')Ar' = G 



Ar' 



\ir')^Ar'.m 



Equations (17) for p{r) and (29) for M{r) allow the 
direct calculation of crs{r) from r = to r = oo. Afterwards, 
we can use Equation (37) to get, via our model, the local 
mean velocity dispersion, ap{R), in the observation plane. 

We have computed both o-g{r) from Equation (40) 
and a-p{R) from Equation (37) for each value of n (n = 
1,2, ...,10) in the Sersic profile (see Fig. 4). In all the 
cases we have calculated the corresponding quadratic mean 
over the total space, i.e. a^ and o-p. We have found that 
{GL/Re)w^ = 3af and a^ = ap with a precision greater 
that the 0.01%, i.e. they are satisfied to at least with four 
significant number. This has been the precision that we have 
used in the calculation of the density, p(r) . 

As we now have a "theoretical" distribution of ap{R), 
we can compute the mean quadratic value, ap{Rm), inside 
any total radius Rm and so we can evaluate the difference 
between the total and any partial integration. In Fig. 5 we 
show the effects that can appear as a consequence of having 



a lack of data in the radial observational interval. We have 
calculated the quadratic mean partial value. 



(^piRm) 



Xf - ap{R)dL{R) 



(41) 



Obviously, when Rm -^ oo we find the total value ap 



1/3(GL/Re 



But when Rm/Re is too low we can find 



important differences between cTp and (jp{Rm), at least for 
some values of n in the Sersic profile (see Fig. 5) . 



5 CONCLUSIONS 

Once we admit that the distribution of the intensity over 
a galactocentric radius of an elliptical galaxy can be well 
represented by the 7?^' " Sersic profile a semi-analytical ex- 
pression for the corresponding spatial density of luminous 
sources, phir), is given. It takes the form of a sum of expo- 
nentials with the same argument, i?^'" ,as in the observed 
intensity profile, I{R)- But for phij-) in each one of these ex- 
ponentials, the argument r^'" is multiplied by a numerical 
factor, Aj, and is easily obtained. Likewise, the correspond- 
ing coefficient of each exponential is easily computed. The 
number {Nap) of these exponentials, the order of approxima- 
tion, depends on the required precision. A number between 
5 and 10 can be sufficient for practical applications. Once 
this semi-analytical expression for the spatial density, p{r), 
is found, the distribution of mass, M{r), the potential, (j>{r), 
and the velocity dispersions, a'1{r) and a^{R), can be com- 
puted in a straightforward manner. 

Furthermore, we show that the total quadratic mean 
of the measured velocity dispersion over the observational 
plane, {GL/R)ap{Rm), cannot take the adequate value 
{GL/R)ap as a consequence of an incomplete integration be- 
cause of the lack of observations. The ratio between ap{Rm) 
and ap obtained by means of computations with parame- 
ters and the function of the "Sersic" models provides us 
with the corresponding correction factor and consequently 
the correct value to use in the virial theorem to deduce the 
mass-to-luminosity ratio. 
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Figure Captions: 



Figure 1: Top: the surface distribution of intensity, 
I{R); Sersic law for diflerent values of the exponent n {n — 
1,2, ..., 10) in Equation (7). The dotted line correspond to 
n = 1. Bottom: the corresponding normalized luminosity 
distribution, L{R)/R, equation (8) and (9). 

Figure 2: Top: plot of the density of luminosity p{r) 
corresponding to the I{R) Sersic profiles for n = 1, 2, ..., 10. 
Bottom: plot of the normalized mass M{r)/M for the same 
values of n. 

Figure 3: Curves for the normalized potential 
(j>{s)/(j){0) for ri = 1,2, ...,10. 

Figure 4: Top: plot of the spatial velocity dispersion 
o'^{r) derived from equation (40) for n = 1, 2, ..., 10. Bottom,: 
plot of the observed velocity dispersion ap{R) derived from 
equation (37) for the same values of n. 

Figure 5: Plot of the aperture velocity dispersion 
o'p{Rm) derived from equation (41) for different values of 
the exponent n (n = 1, 2, ..., 10) in the Sersic profile. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science MJfjX style file. 



© 0000 RAS, MNRAS 000, 000-000 







■5 



K 



oo -10 

o 



15 



20 



0.8 



^ 0.6 
^ 0.4 



0.2 





4 



2 2 4 

log (R/Re) 



o 







■5 



10 



15 



20 




4 



2 2 

log (r/Re) 



4 




4 



2 2 4 

log (r/Re) 



"1 I r 



1 I I I I I r 



0.8 



o 



0.6 



^ 0.4 



0.2 








n=10 



J I I I I I L 



J I L 



4 



2 2 

log (r/Re) 



4 



0.8 



0.6 



b 0.4 



0.2 








J I I I I I I I I I L 



4 



-2 

log (r/Re) 



2 



K 



0.8 



0.6 



b 0.4 



0.2 








J I I I I I I I I I L 



4 



-2 

log (R/Re) 



2 



"1 I r 



"1 I r 



"1 I r 



0.8 



0.6 



K 



b 



0.4 



0.2 







J I L 



4 



n=10 




J I L 



J I L 



2 

log (R /Re) 



2 



